Magnetization damping in a local-density approximation 
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The linear response of itinerant transition metal ferromagnets to transverse magnetic fields is 
studied in a self-consistent adiabatic local-density approximation. The susceptibility is calculated 
from a microscopic Hamiltonian, including spin-conserving impurities, impurity induced spin-orbit 
interaction and magnetic impurities using the Keldysh formalism. The Gilbert damping constant 
in the Landau-Lifshitz-Gilbert equation is identified, parametrized by an effective transverse spin 
dephasing rate, and is found to be inversely proportional to the exchange splitting. Our result 
justify the phenomenological treatment of transverse spin dephasing in the study of current-induced 
magnetization dynamics in weak, itinerant ferromagnets by Tserkovnyak et al.} We show that 
neglect of gradient corrections in the quasiclassical transport equations leads to incorrect results 
when the exchange potential becomes of the order of the Fermi energy. 

PACS numbers: 72.25.Rb,75.45.+j,76.50.-l-g 



I. INTRODUCTION 

The drive to miniaturize and reduce power demands 
of electronic appliances motivates research in nanoscale 
magnetoelectronics, i.e. the science and technology that 
exploits additional functionalities offered by ferromag- 
nets integrated into electronic circuits and devices. Spec- 
tacular advances have been realized already in the last 
decade, mainly in the area of magnetic disk and magnetic 
random access memories. Itinerant transition metals and 
its alloys are the materials of choice for magnetoelectronic 
applications due to their high electric conductivity and 
Curie temperatures. Increasing speed and reducing en- 
ergy demands of switching a bit of information encoded 
by the magnetization direction of a ferromagnetic grain 
is one of the key problems in the field. A thorough under- 
standing of the dynamics of the magnetization order pa- 
rameter in transition metals is necessary to make progress 
in this direction. 

Phenomenologically, the low temperature magnetiza- 
tion dynamics in ferromagnets is well described by the 
Landau-Lifshitz-Gilbert (LLG) equation. ^'^ Ferromag- 
netic resonance (FMR) experiments can be fitted to ob- 
tain accurate values for the parameters of the LLG equa- 
tion, viz. the Gilbert constant that parametrizes vis- 
cous damping and the effective (demagnetization and 
crystal anisotropy) fields. The LLG phenomenology has 
been successfully applied to explain a rich variety of dy- 
namic magnetic phenomena.'^ Some progress has been 
made in predicting magnetic crystal anisotropics by first- 
principles calculations.® However, in spite of being a cru- 
cial device parameter that governs the switching time of 
magnetic memory elements, the material dependence of 
the intrinsic magnetization damping has not yet been un- 
derstood. The Gilbert damping parameter also plays an 
important role in current-induced magnetization excita- 
tions and domain wall motionfiii 



Deriving a microscopic description of the dynamics of 
transition metal ferromagnets is a formidable task; even 
the nature of the ground state is still under debate. Two 
different viewpoints can be distinguished. On one hand, 
ferromagnetism can be seen to be caused by the atomic 
correlations in partially filled and essentially localized d 
orbitals. The s-electrons that are responsible for electron 
transport arc in this picture affected by the magnetic or- 
der only indirectly via local exchange interactions. Such 
physics is expressed by the so-called s-d model in which 
the spin of localized c?-electrons. Si, interact with free 
s-electron spins, s(ri,t), through a Heisenberg exchange 
term. 

In the opposite point of view, the d-electrons are not 
only broadened into bands, but are also strongly hy- 
bridized with neighboring s-p orbitals. A separate treat- 
ment of the orbital symmetry is then not warranted for 
the description of low energy properties at long time 
scales. The Stoner model represents the essence of this 
itinerant magnetism in terms of two (minority and ma- 
jority) parabolic energy bands that are split by a con- 
stant exchange potential. Spin density-functional theory 
in a local density approximation is the modern version 
of itinerant magnetism, forming the basis of most band 
structure calculations to date. The nature of the real 
wave function of 3d ferromagnets that combines features 
of both extremes is presumably captured by sophisticated 
many-body frameworks such as the dynamical mean-field 
model. It is at present not obvious, however, how to com- 
pute the low-energy collective dynamics of ferromagnets 
taking disorder as a well as local correlations into ac- 
count. 

The Gilbert damping coefficient in the LLG equation, 
usually denoted by a, has attracted quite some theoreti- 
cal attention. Incoherent scattering of electron-hole pair 
excitations by phonons and magnons is a possible mech- 
anism by which energy and angular momentum can be 
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dissipated. Heinrich et alM- suggested a model in which 
conduction electron spins become polarized by scatter- 
ing with magnons. The spin angular momentum is sub- 
sequently transferred to the lattice by spin-orbit medi- 
ated relaxations. The resulting damping coefficient was 
found to be proportional to the electronic scattering rate, 
a ~ T^^. We will return to this result in Sec. lIVBl More 
recently, a phenomenological treatment of the Gilbert 
damping has also been reported in Ref. H. 

A different relaxation process was proposed in Ref. ITol . 
and further elaborated in Rcfs. 11,12,13.14: In the pres- 
ence of spin-orbit interaction, the electronic energy lev- 
els depend on a time-dependent magnetization direction, 
giving rise to the notion of a "breathing Fermi surface". 
The time lag of the electronic distribution response to a 
moving magnetization vector is equivalent with dissipa- 
tion. In this model the Gilbert damping coefficient is pro- 
portional to the scattering time, a ^ t. Extrinsic contri- 
butions to the FMR linewidth such as eddy currents ex- 
cited by time-dependent magnetic fieldsfi^ sample inho- 
mogeneities, or two-magnon scattering processe o^^'^^i^^ 
have been suggested as well. 

In diluted magnetic semiconductors such as 
(Ga,Mn)As the magnetism originates mainly from 
the local spins of the half-filled spin-5/2 Mn d-shells. 
The spins are coupled by a local exchange interaction to 
the valence band holes, and non-locally, via the holes, 
ferromagnetically to each other. The holes contribute 
only slightly to the magnetization, but are exclusively 
responsible for the finite conductivity. The s-d model 
is therefore appropriate for understanding the magne- 
tization damping in ferromagnetic semiconductors f^^*^ 
Magnetization damping in the s-d model can be 
understood in terms of the so-called spin-pumping 
mechanismi2ii22ii2^i2i The motion of the localized spins 
pumps a spin current into the conduction electron 
bath, in which the thus created spin accumulation is 
dissipated by spin-flip scattering. Ref. HO reported a 
non-monotonous dependence of the damping on the 
scattering rate, i.e. a ~ for clean and a ^ t for dirty 
samples. As mentioned above, the s-d model does not 
necessarily give a good description of transport and dy- 
namical properties of transition metal ferromagnets. The 
notion of d-electrons pumping spins into an s-electron 
system becomes doubtful when the hybridization is very 
strong. Recently it has been demonstrated that the 
magnetization dynamics in the s-d model and in an 
itinerant Stoner model can be quite different indeedji 
For example, for a given spin-flip relaxation mechanism, 
the Gilbert damping is significantly suppressed in the 
s-d description by a factor of the (usually small) fraction 
of the total magnetization carried by the delocalized 
s-electrons. 

Recently Kohno et al^ reported results for both the 
Gilbert damping and the non-adiabatic current-induced 
spin torque term postulated by Zhang and Lee^ which is 
parameterized by a constant /3. The results in Ref. 1 have 
been obtained under the assumption that the exchange 



splitting is small compared to the Fermi energy. The dia- 
grammatic treatment by Kohno et al. was not restricted 
to weak ferromagnets, but the differences turned out to 
be very small for transition metals.- 

In the present paper we generalize the treatment of the 
transverse spin dephasing of Ref. U beyond the relaxation 
time approximation. We relax the previous limitation to 
weak ferromagnets and derive the corresponding Gilbert 
damping. We use a self-consistent adiabatic local-density 
approximation (ALDA) model in the presence of a di- 
lute concentration of scalar and magnetic impurities, as 
well as spin-orbit interaction originating from impurities, 
and we demonstrate how to generalize the previous treat- 
ment to strong ferromagnets. The generalization of the 
Keldysh approach is non-trivial and introduces subtle, 
but important gradient corrections. 

Our main result is that for spatially homogeneous 
itinerant ferromagnets, the Gilbert damping constant is 
given by 



where A is the modulus of the local-density exchange- 
correlation potential, and tJ^^ is a transverse (Bloch) 
spin dephasing rate caused by spin-orbit interaction and 
magnetic disorder. This appears to be at variance with 
Kohno et al.,^^ who find a Gilbert damping constant that 
depends on both the longitudinal and the transverse scat- 
tering rates. Except for this issue, we obtain the same 
detailed expression for a. This is gratifying since these 
theoretical machineries are completely different. 

This paper is organized in the following way: The mi- 
croscopic model, as well as the simplifying ALDA are 
presented in Sec. [Til while the linear response formalism, 
microscopically and phenomenologically, will be treated 
in Sec. IIIII The detailed derivation of the linear response 
function, starting from the Keldysh Green function for- 
malism, is the topic of Sec. IIVI Conclusions are summa- 
rized in Sec. El 

II. TIME-DEPENDENT ADIABATIC 
LOCAL-DENSITY APPROXIMATION (ALDA) 

Density-functional theory (DFT) is a successful and 
widely used method to study electronic structure and 
magnetism in transition metal ferromagnets In the 
Kohn-Sham implementation, non-interacting pseudo- 
particles are introduced that exhibit the same ground 
state density as the interacting many-electron system. 
This is realized by introducing a fictitious exchange- 
correlation potential that has to be determined self- 
consistently by energy minimization?2I DFT can be ex- 
panded to handle time-dependent phenomena in systems 
out of equilibrium.— 

We study the magnetization dynamics in a simplified 
time-dependent DFT in the local spin-density approx- 
imation, in which non-interacting Kohn-Sham particles 
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are treated as free electrons. We realize that the transi- 
tion metals have complex energy bands and wave func- 
tions also in the local density approximation and that 
even integrated properties such as conductivities are not 
necessarily well described by free electrons. However, 
since we are interested in dirty systems, in which an ad- 
ditional scrambling occurs by elastic impurity scattering 
between different bands, we are confident that our treat- 
ment is a good starting point for more sophisticated com- 
putations that take the full band structure into account. 

Our model for itinerant ferromagnetism is described by 
the Hamiltonian 

H = i [Ho + V{r)] + ]pnHa, + ]^-fhH^,[p{r,t)] ■ a 

+ ]pnh{r, t)-a + VUr) + K„(r). (2) 

Matrices in 2 x 2 spin space are denoted by a hat ("). 
Spin- independence is indicated by the unit matrix 1, and 
o- is a vector of the Pauli matrices. Here, Tig is the 
translationally-invariant Hamiltonian of non-interacting 
electrons, and V{r) is the elastic spin-conserving impu- 
rity potential. iJ > is an effective magnetic field in 
the z-direction consisting of internal anisotropy fields and 
externally applied contributions, and —7 < is the elec- 
tronic gyromagnetic ratio. Electron-electron interactions 
are described by the exchange-correlation (vector) field 
ffxc- The weak, transverse magnetic driving field is de- 
noted h{r,t), and the potentials due to impurity induced 
spin-orbit interaction and a magnetic disorder configura- 
tion are denoted Vso{r) and K„(r), respectively. 

The published approximation schemes for time- 
dependent exchange-correlation functionals are still 
rather crude and/or untested. The simplest approxi- 
mation is known as the adiabatic local-density approx- 
imation (ALDA). Here the conventional LDA exchange- 
correlation potential is adopted for the instantaneous 
time-dependent densityi^i^S The ALDA potential is 
therefore local in both spatial and temporal degrees of 
freedom and reduces for the current problem to 



'yhH^^[p]{r,t) « Am{r,t), 



(3) 



where A is an effective exchange splitting constant and 
m{r,t) is the local magnetization direction of the ferro- 
magnet. By construction, the exchange field is always 
parallel to the magnetization direction, and thus auto- 
matically satisfies the zero-torque theorem<^ Since we 
are interested in the low-energy transverse magnetization 
dynamics the (atomic-scale) position dependence of A is 
disregarded. 

In the ALDA, the exchange-correlation potential is re- 
leased from a possible functional dependence on the his- 
tory of the system. The ALDA should therefore be valid 
only when the system is close to the equilibrium configu- 
ration, i.e. for a slowly varying magnetization direction 
in both space and time. This is the case when hdt <C A, 
and dr ^ kp, with kp being a characteristic Fermi wave 



vector. Improved descriptions of the exchange potential 
have been proposed {e.g. the generalized gradient ap- 
proximation), but for slowly varying uniform perturba- 
tions, such corrections are believed to be smalli^ 



III. LINEAR RESPONSE 

For a weak magnetic driving field, the response of the 
ferromagnet can be formulated within linear response 
theory. An expression for the response to the per- 
turbative field h{r,t) is derived quantum mechanically 
from the ALDA Hamiltonian defined by Eqs. ^ [3]) in 
Sec. nil Al The response derived from the phenomcno- 
logical Landau-Lifshitz-Gilbert equation is presented in 
Sec. IIIIBI These results are then used in Sec. |IV]to find 
a microscopic expression for the Gilbert damping coeffi- 
cient. 



A. Quantum linear response 

The Kubo formalism provides expressions for the lin- 
ear response to a time-dependent perturbation. The re- 
sponse functions can be derived by considering the time 
evolution of the non-equilibrium density matrix. Start- 
ing from the effective Hamiltonian ^ in the ALDA of 
Eq. ([3]), the small time-dependent perturbation operator 
should include the self-consistent exchange as: 



Hintit) ^ Jdr 



■-fh{r,t) 



s(r), 



where s{r) is the spin density operator. The empha- 
sis of this paper is on the transverse, non-equilibrium 
components of the spin density. We denote sq = \sq\, 
(s) = — sqCz + (<5s), where (Ss) ± e^. Hence, \s\ = sq, 
and m — —{s)/so in the ALDA. 

For axially symmetric systems, the non-equilibrium 
spin density response can be expressed conveniently in 
terms of 6s± — dsx ± iSsy. The transverse part of the 
response to the magnetic field can then be written as 



A 



(4) 



where the retarded susceptibility tensor 

X^Ar,r';t) = ^e(i)([s^(r, t), s.(r', 0)]), 

has been introduced. The brackets [. . . ] indicate a com- 
mutator, and the angular brackets (...) a thermody- 
namical average. In the derivation of the above expres- 
sion, we have made use of axial symmetry under which 
X++ = — 0. In the ALDA, Eq. ^ can be simphfied 
to 



{Ss^{q,u)) = -x-+(g,tj)/i_(q,w). 



(5) 
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where the self-consistent hnear response to the driving 
field 



which corresponds to a susceptibility 



A 



(6) 



has been introduced. Hence, in the ALDA, the linear 
response of an interacting system reduces itself to calcu- 
lating the response x h of ^ noninteracting system with 

a fixed (Stoner enhancement) exchange fieldi^ 



B. Landau-Lifshitz-Gilbert susceptibility 

The phenomenological Landau-Lifshitz equation^, is 
widely used to model transverse magnetization dynam- 
ics. The magnetization direction m{r,t) is treated as a 
classical field whose dynamics is governed by an effective 
magnetic field Hcs{r, t), obtainable from the free-energy 
functional of the system, F[iVf]: 

H,s{r,t)^-dMF[M]. 

The Landau-Lifshitz equation describes undamped [i.e., 
free-energy conserving) precessional motion about the lo- 
cal effective magnetic field: 

dtm{r,t) = -■ym{r,t) x Hcs{r,t), 

preserving the magnitude of the magnetization. The field 
Hc.f[{r, t) includes contributions from external, exchange, 
demagnetization and crystal-anisotropy magnetic fields. 

The Landau-Lifshitz equation does not dissipate en- 
ergy since the effective magnetic field always points nor- 
mal to the instantaneous constant-free-energy surfaces. 
However, the electronic degrees of freedom do not re- 
spond infinitely fast to the magnetization dynamics, 
which means that in reality the effective field is a func- 
tional of the time-dependent magnetization at previous 
times. A finite lag in the response of the dynamics corre- 
sponds to energy dissipation. In a magnetic system, the 
energy-loss implies a lowering of the Zeeman energy by 
a torque in the direction of the cross product of magne- 
tization and its time derivative; the energy-loss can be 
parametrized by the phenomenological Gilbert damping 
constant Hence, we arrive at the Landau-Lifshitz- 
Gilbert (LLG) equation: 

dtiri = — 7Tn x Hcs + am x dtm. 

Here, Hcs only depends on the instantaneous magnetic 
configuration of the ferromagnet. Generally, the damping 
is a tensor quantity with symmetries reflecting the crys- 
tal structure, but in practice anisotropic corrections to 
damping are small compared to those in the free-energy.'^'* 
Assuming an external field of the form jHcsir.t) = 
ujo{r,t)ez, and a small rf driving field h{r,t), the ex- 
cited small-angle transverse magnetization dynamics can 
be computed easily by the linearized LLG equation. 



m_(q,w) 



7/i_(q,w) 



ISO 



uJo{q, Lu) ~ Lu ~ ia(<7, lu)uj ' 



that can be directly compared with the microscopic re- 
sponse function x h by Eq. ([6]). Assuming that A//i ^ 

{u!,u!q) one obtainsiSiSS 



a{q, uj ^ 0) = 



A2 



7^2 So 



lim d^lmx-+{q,uj) 



(7) 



Hence, finding a microscopic expression for the Gilbert 
damping is equivalent to determining the quantum me- 
chanical transverse susceptibility. 

It is worth noticing that, in general, the damping co- 
efficient may depend on the spin-wave wave vector q. 
A damping of the form a(q,w ^ 0) ~ will intro- 
duce an additional dissipativc term in the LLG equa- 
tion, dtm oc —am x V^dtm, which is similar in form to 
the exchange field oc m x V^m for anisotropic fcrromag- 
nets. The main emphasis in this paper, however, is on 
the isotropic part of the Gilbert damping. 



IV. MICROSCOPIC DERIVATION OF THE 
SUSCEPTIBILITY 

In this section we determine the susceptibility function 
according to time-dependent spin DFT in the ALDA for 
a disordered ferromagnet. The Keldysh Green function 
formalism is used, and the assumption of weak and slowly 
varying perturbations in space and time allows us to use 
the simplifying gradient expansion. Finally, Eq. ([7| is 
invoked to obtain the Gilbert damping coefficient. 



A. Kinetic equation 

We proceed from the ALDA Hamiltonian (see Sec. HI)): 

H^l[Ho + V{r)] + ^(A + jhH)a, 

+ ^-rhhir,t)-a + V,,{r) + Vn.ir). (8) 

In the following discussion, we assume a homogeneous 
static magnetic field H, define A' = A -I- jhH, and drop 
the prime for brevity. 

The impurities are assumed to be randomly distributed 
over positions with short-range, scalar disorder poten- 
tials 



a;o(<7, uj) — Lo — ia{q, U!)lu 



The scattering potentials are Gaussian-distributed with 
zero average and a white noise correlator 

{V{r)Vir')) =^S{r~r'). 
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We define a characteristic scattering time r by = 
7r(i/| + v^)T/h, with Vg being the density of states at the 
Fermi level for electrons with spin s. The spin-orbit in- 
teraction associated to impurities is described by the po- 
tential 

%o[r)^iP&-{VV{r)xV), 

where /9 is a spin-orbit interaction strength given by 
—h^/Amlc^, in terms of the electron mass nie and the 
speed of light c. The magnetic disorder in the ferromag- 
net is modeled as 

Vmir) = ^Vm{ri)S{r - ri)S{ri) ■ a, 



where S{ri) denotes the spin of an impurity at position 
r.i. The internal degrees of freedom of the magnetic im- 
purities are assumed to be frozen. Also the vector impu- 
rity exchange field V(t') = iTr[V^n ('")''■] is taken to be 
distributed according to Gaussian white noise character- 
istics, i.e. {Va{r)) = and 

where a, (3 denote spatial components of the vector field, 
and the strength of the second moment is given by 



Sm 



x,y 

2, 



similar to Ref. |25!. 

We employ the Keldysh Green function formalism^^ 
to calculate the spin susceptibility defined earlier. This 
method has distinct advantages over the equilibrium for- 
malism when it comes to describing non-equilibrium phe- 
nomena, but also some drawbacks. However, as we will 
see in Sec. IIVBI in particular, whereas the equilibrium 
formalism requires tedious calculations of vertex correc- 
tions, the Keldysh formalism requires that one carefully 
accounts for subtle gradient corrections in order to obtain 
the correct dynamics. 

The Green function in Keldysh space, denoted by an 
inverted caret ( " ) takes on the form'^^ 



^, , /G«(l,2) G^(l,2) 
^^^'^^~(v G^(l,2) 



The retarded, advanced and Keldysh Green functions are 
given by 

G«(l,2) = -ie(ti-t2)({*(l),*n2)}), 
G^(l,2) = +19(^2 -ti)({*(l),*n2)}), 



and 



G^(l,2) = -i([*(l),vl/t(2)]), 



respectively. Brackets [. . . ] indicate a commutator, the 
curly brackets {••.} an anti-commutator, while ^''^^^ is 



a fermion annihilation (creation) operator. In this no- 
tation, all field variables (position, time and spin) are 
contained in the numerical indices 1 and 2. 

In the presence of slowly varying perturbations, the 
two-point propagator variables can be transformed into 
the Wigner representation, viz. the center-of-mass coor- 
dinates and the Fourier transform of the Green function 
with respect to the relative coordinates: 



G(X, k) 



Ax e 



'ikx 



G{X + x/2,X -x/2). 



Here, a four-vector formulation has been introduced, 
where the vector for the center-of-mass coordinates is 
X = {R, T) , the corresponding relative coordinates are 
given by X = (t',^), and finally k — {k,e). The four- 
vector product is defined —et+k-r. The Wigner 
representation is particularly convenient when the varia- 
tion of the Green function on center coordinates is slow 
on the scale of the Fermi wavelength, since this allows 
us to perform a gradient expansion in these coordinates. 
Subtracting the Dyson equation and its conjugate, one 
finds the relation 



[Go-i - S ? G] = 0, 



(9) 



where the symbol ® denotes a convolution (in position, 
time and spin), the commutator corresponds to the 2 x 
2 Keldysh matrix structure, and Gq^ is the inverse of 
the retarded/advanced Green function in absence of any 
impurities. It is diagonal in Keldysh space, and in the 
Wigner representation has the structure 



^Aa,-^-fhh{R,T)-& 
(10) 

in spin space, with e^, denoting the free-electron energy 
measured with respect to the chemical potential. The fi- 
nal component in Eq. ([9]), E, is the self-energy due to the 
impurity configurations and spin-orbit interaction. One 
can show by a formal Taylor expansion that the convo- 
lution can be represented by 

{A (g, B){X, k) = e^'^^^-^^-^^-^^'^/^AiX, k)B{X, k) (11) 

in the Wigner representation!^ 

Physical quantities such as occupation probabilities 
and densities are expressible in terms of the distribution 
Green function G^, which is given by the combination 



G^ 



1 



G 



K 



where we have introduced the spectral function A — 
\{G^ — G^). To derive a kinetic equation for G^, we 
subtract the diagonal components of Eq. ([5]), and com- 
bine the result with the Keldysh component of the same 
equation. In summary, one finds the kinetic equation 



G' 



>G< -G^ 



G' 



= S< (g) G^ - G' 
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Assuming slowly varying perturbations, we now use the 
gradient expansion, in which the exponential in Eq. (jlip 
is expanded, and only the first two terms of the expansion 
are kept.^^ This results in a simplified kinetic equation 
for the distribution Green function, viz. 

= A([S«,G%-[G<,S-4]p) 

-i([G^S<]p-p<,G^]p). (12) 

All terms to first order in the generalized Poisson bracket 
are kept, [X, Y]p = dxX ■ dkY — dkX ■ dxY, where the 
four- vector notation implies dx • dk = Or • dk — hdrd^. 
We see that the gradient expansion reduces the complex 
convolution of the Dyson equation ([9]) to the 2x2 matrix 
multiplication of Eq. (|12p . 

It is important to carefully include all contributions 
to first order in the Poisson brackets. The correct dy- 
namics is captured only when gradient corrections to the 
spectral function are kept, and, in the case of anisotrop- 
ically distributed magnetic impurities, those to the self- 
energies as well. Such gradient corrections are caused by 
the non-uniform driving field. For simple systems, e.g. 
normal metals, such corrections are absent. This means 
that one should not expect naive generalizations of Boltz- 
mann equations to strong ferromagnets to be completely 
correct. 

The final ingredient that transforms Eq. into a 
useful kinetic equation is an expression for the self-energy. 
For weak impurity scattering, the self-consistent Born 
approximation is appropriate: 

E(l,2) = (ftot(l)G(l,2)ytot(2)), 

where Vtot is the total potential (which is diagonal in 
Keldysh space), angular brackets (...) denotes impurity 
potential averaging, while G(l,2) is already an impurity 
averaged Green function. This expression can be sepa- 
rated into four different self-energy contributions. The 
spin-conserving impurity scattering is described by 

/dk' 
^G{R,T-k',e). 

Introducing n = k x k' , the terms arising from the spin- 



orbit interaction can be written as 

E« (H, T; fc, e) = -i^P J ^ [G(fi, T; k' , e)& ■ n 

-n-&G{R,T;k\e)], 

and 

/c]k' 
j^n-aGiR,T;k',e)&-n. 

The magnetic impurity configuration results in 
S„,(i2,r;fe,e)= J2 J ^^rGiR,T;k',e)a,. 

i—x^y^z 

Finally, to make connection with the spin density, one 
can use that 

All necessary quantities are now defined, and a kinetic 
equation for the distribution Green function can be de- 
rived. In the next section, the details are worked out for 
a bulk, single-domain ferromagnet. 



B. Homogeneous ferromagnet 

We concentrate in the following on the Gilbert damp- 
ing constant in the limit of vanishing spin-wave wave vec- 
tor, q — > as measured in FMR experiments. In this 
limit (and the ALDA), only spin-orbit interaction and 
magnetic disorder can transfer angular momentum out 
of the spin dynamics into the lattice. Without it, spin 
and orbital degrees of freedom are completely decoupled 
and the Gilbert constant vanishes. Spin waves with finite 
wave lengths may decay also by spin-conserving scatter- 
ing, which is likely to dominate magnetic impurity or 
spin-orbit interaction scattering when q becomes larger 

We simplify the notation by defining the time and 
energy-dependent density matrix 

/r\k - 
j^G<iT;k,e), 

and solve Eq. p2)l to obtain a diffusion equation for this 
quantity. In detail, we find that 
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+ E^^'/ / tS^"'"^ (a.G%G< -G«a.G<a, - (G<a.G^a, - a.G<a,G^)) 
+ /(^/ (|^«»".([-^G%,G<],-[G«,a.G<a,]p-([G<,a^ 



Here, the arguments to p{T,e, G^^^{T;k,£) and 
G^{T; k' , e) are not written out explicitly for the sake 
of notation. Summation indices i,j run over Cartesian 
components x, y, z. On the left hand side of the equa- 
tion, we recognize precession around the fixed exchange 
field and the driving field, as well as a gradient term due 
to the non-uniformity of the driving field. On the right 
hand side we find collision integrals due to spin-orbit in- 
teraction and magnetic impurities. We also find gradient 
corrections to the collision integrals in the third and fifth 
line above. These corrections are often neglected, but are 
important for strong ferromagnets, to be discussed below. 
As explained, scalar disorder does not affect the uniform 



spin dynamics and drop out of the kinetic equation. 

For the response function x h(^) introduced in 

Sec. IIII Al we need to find an expression for {Ss-{uj)), 
the transverse part of the spin density. To this end, we 
extract the upper right matrix component of p(T, e) , a 
matrix component we simply denote ^s„(T, e). This is 
now related to the non-equilibrium spin density by 



(5._(T)) = 




de 5s_(r,e). 



With i^{ek) being the density of states at energy Sk, we 
find that 



ihdTSs-{T,e) - ASs-{T,e)+jhh-{T)p4e) + -jtfdTh^iT) J det i^{ek)d,G<{k,e) 
= - 2i(e± +^||) Jdsk Hek) [A4k,e)ds^iT,e) ~ p4e)A^iT;k,e)] 

-HU-^\\) J dsk Hek)[deA,{k,e)dTSs.{T,e) - d,p,.{e)dTA.{T;k,e)] 

-4(e±-e||) J dsfeKefc) [ReGf{k,s)Ss.{T,e) + p,{e)ReG^{T;k,s)] (14) 
+ 2i?i(a+^||) J dekiyiek) [d,pd{e)dTReG'^iT;k,e) + deReG^{k,e)dTSs-iT,e)] 
- 1^/3^ I dsk Hek)k^ I d£fc' 'y{ek')k'^ [Ad{k, e)G<{T; k', e) - G< (fc', e)A^(T- k, e)] 
+ ^^I3^ Jdek i^{ek)k^ J dsk' y{ek')k'^ [9,G<(fc', £)5TReG«(T; ft, e) + 9,ReGf (ft, £)9tG<(T; fc', e)] , 
where have used the convenient matrix notation 

G^lGd + G&. (15) 

We now need to Fourier transform and integrate this formidable equation over energy to obtain a diffusion equation 
for (5s_(w)). Before we proceed with this calculation, we notice that the real part of the response function simply 
determines the resonance condition for the system and is thus unimportant for the determination of the Gilbert 
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damping. With this in mind, we write 



(^iw- A)(,5s_(tj)) = -(e_L+C||) / de / dsk iy{ek)[pd{£)A-{uj;k,e) - Ad{k,e)6s_{uj,e)] 



(16) 



/ de / d£feK£fc)fc' / det' >^{ek')k'^ [G<{k' ,e)A^{Lj;k,s) ~ A4k,e)G<{u;;k\e)] +T{Lo)h^{u), 



r 



where J-{uj) is real, and thus does not contribute to 
the damping. The dissipation is now determined by 
the above integrals over energy and momentum. No- 
tice how the signs between the longitudinal and trans- 
verse magnetic impurity scattering strength enter in the 
above equation. For simple isotropic magnetic impuri- 
ties, i.e. with = , the second line of Eq. (fTB|) does 
not contribute to the dissipative dynamics. This term is 
due to gradient corrections involving self-energies in the 
original kinetic equation (|12p . To correctly capture the 
dynamics when the magnetic impurities are anisotropi- 
cally distributed, it is essential to include such gradient 
corrections as well. 

In order to calculate the integrals, we need expressions 
for , (5s_ and the spectral function A. Since these 
quantities enter the collision integrals, we can solve for 
5s- and A to zeroth order in scattering rates. To this 
end, we solve Eq. for the Fourier transform of , 
and find that 

G<i.-,k'^e)^^-^(l^'^)G<i^k',e) 



2A 



a,G<(fc',e)-HO(e,e<:'). (17) 



At this point we take G^{k',e) — mF{£)Ad{k',e), with 
np the Fermi-Dirac distribution, so that 

d,G< (fc', e) - -i(5(£)Ad(fc', e) + mF{£)dMk', e) 

at low temperatures. Additionally, we use that 

dk' 



Ss- (w, e) 



:G<{oj-k',e), 



in combination with Eq. (jl7p to solve for the second quan- 
tity in question. 

In the dilute limit, the Lorentzian shape of the spec- 
tral function approaches a Dirac delta function, and two 
quasiparticle spin bands, split by the exchange field A, 
are resolved. For a uniform, time-independent transverse 
magnetic field one finds 

io(fc,e) =7r^5(e-efe,) M + sg, + s ^^^' j . (18) 

The two spin bands are denoted by s =T, i= ±, and the 
notation Eks = £fc + sA/2 has been introduced. A non- 
uniform driving will also introduce terms in the spectral 



function that are linear in gradients. A detailed deriva- 
tion of the spectral function to first order in Poisson 
brackets for a time-dependent driving field h[T) is given 
in App. 1^ with the main result being 

A{T-k,e) = Ao{T;k,e) 

+ ^^{^Mk,£) + ^dMk,e)^a,dT{h-a), (19) 

where Aq is the spectral function from Eq. (HH]). We see 
that the weak, transverse driving field induce off-diagonal 
gradient corrections to the "instantaneous" spectral func- 
tion. The diagonal components are unchanged, and are 
given by 



and 



Ad{k,e) = 7r^(5(e - £ks 



Az{k,£) =TT^s5{e-eks)- 



We are now in a position to calculate the above energy 
integrals. To be more specific, considering the integrals 
due to magnetic disorder, we find that these terms be- 
come 



defc v{£k) [pd{£)A^{uj]k,e) - Ad{k,e)5s^{uj,£)] 



2A 



-(5(e) (i^T + n) 



and 



ftw J dekviek)[deA2{k,e)Ss^{uj,e) - dePz{£)A^{uj:k,e) 

= ^ 5{£){v^ -vyf. 

One can derive analogous results also for the spin-orbit 
contribution in Eq. (jl6p : inserted we find that 



[hw - A)((5s_(w)) « T{uj)h^ 
2iTT'yfr' 
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Using Eq. Q to identify (the low frequency) Imx |_, w( 

find from Eq. ([7]) that the Gilbert damping is given by 



2'Kh 
9^ 



So 



[Cil^^T' + ^D + S^II'^Tn] , (20) 



which agrees with the diagrammatic calculation of 
Kohno et al.^ 

Next, we would like to relate the Gilbert damping con- 
stant in Eq. (j20p to other physical quantities. Comparing 
a Bloch-Bloembergenii^S equation of motion for the mag- 
netization vector with the corresponding LLG equation, 
we find for weak driving fields and small angle magneti- 
zation dynamics that 



Atj_' 



(21) 



where, in our case. 



Above we have defined the effective transverse scattering 
rates from spin-orbit interaction and magnetic impuri- 
ties, viz. 

and 

;^ = ^[a(^T+n') + 2e|inn]- 

By comparison, the longitudinal spin relaxation rate ob- 
tained from e.g. Fermi's Golden rule reads 



1 47r^ 



For weak ferromagnets, the density of states and momen- 
tum at the Fermi energy is not strongly spin dependent, 
i.e. Vs ~ vp, and kps — kp. Therefore 2so ~ HAvp, 
which implies equal transverse and longitudinal scatter- 
ing rates for impurity induced spin-orbit interaction and 
isotropic magnetic impurity scattering, i.e. t± = t\\.^^ 

Since we succeeded in reproducing the general dia- 
grammatic result of Ref . we also identified the neces- 
sary measure to transcend the semiclassical treatment of 
Ref. m. Most important are the gradient corrections to 
the spectral function, but in the presence of anisotrop- 
ically distributed magnetic impurities, gradient correc- 
tions to the self-energies should be included as well. 

In a metal, the longitudinal spin-orbit induced scat- 
tering time depends on the spin-conserving elastic scat- 
tering time. Experimentally, one typically finds that the 
ratio of spin-conserving to non-spin-conserving scattering 
events e = t/tm is not very sensitive to the concentration 



of impurities. ^''^ This means that in systems where spin- 
orbit induced dephasing dominates the Gilbert damping, 
a is proportional to the resistivity of the system. 



hne 
2Am 



where n is the electron concentration, and p is the re- 
sistivity. A linear relation between Gilbert damping and 
resistivity was found in a recent experimental study of 
electronic transport in thin permalloy films by Ingvars- 
son et al^ 

The present discussion is mainly focused on the low 
temperature regime, but the linear dependence on the 
scattering rates suggests that the damping should in- 
crease with increasing temperatures. It is experimen- 
tally known that e is not very sensitive to variations in 
temperaturci^ Consequently, we expect that the Gilbert 
damping constant is proportional to the resistivity also 
at higher temperatures. 

Since the Gilbert damping coefhcient is proportional 
to A~^, dissipation is reduced in the strong ferromagnet 
limit. Since the damping in Eq. (|21[) depends on the 
transverse spin dephasing rate, Gilbert damping does not 
vanish in half metals in which the chemical potential falls 
below the band edge of one of the spin bands. 

The form of the damping coefficient found in 
Eqs. (POI [^T|l agrees with previous studies on spin-flip 
mechanisms for magnetization damping. Nearly four 
decades ago, Heinrich et al^ suggested, based on the s- 
d exchange interaction between localized d electrons and 
itinerant s electrons, that electron-hole pairs could be ex- 
cited by magnons. Assuming that the exchange splitting 
is much larger than the spin-flip rate, and denoting the 
fraction of the total spin carried by the delocalized elec- 
trons by 77 < 1, they found a = rjh^vp /2sQTsf in the long 
wave length limit. The result is expressed in terms of 
Tsf , a phenomenological electron-hole pair life time. This 
result can be compared with Eqs. (^01 [^Tjl by using the 
approximation 2sq « hAup. We see that when the mag- 
netization is mainly carried by the d-electrons, which are 
not affected by spin-flip scattering, the predicted damp- 
ing in the s-d model is much weaker than in our Stoner 
model. 



V. CONCLUSION 

In conclusion, we present a kinetic equation for the 
distribution matrix of itinerant ferromagnets in the adi- 
abatic local-density approximation. The spin suscepti- 
bility and Gilbert damping constant is obtained micro- 
scopically for a homogeneous ferromagnet by the Keldysh 
Green function formalism. Magnetization damping arises 
from magnetic disorder in the ferromagnet, and we have 
shown that it is important to keep all terms to linear 
order in Poisson brackets to obtain the correct result in 
presence of impurity induced spin-orbit interaction mag- 
netic disorder. The Gilbert coefficient can be expressed 
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in terms of an effective transverse spin dephasing rate 
that has been introduced earher as a phenomenologi- 
cal constantji. Our framework could also be generalized 
to handle first-principles band structure calculations for 
specific types of impurities and disorder, that may be 
the starting point for systematic studies of the Gilbert 
damping as a function of material parameters. 



with det(. . . ) denoting a matrix determinant, is simply 
the inverted retarded Green function to zeroth order in 
Poisson brackets. Matrix manipulations results in 



G«(T;fc,e)«G^(r;fc,£) 
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APPENDIX A: SPECTRAL FUNCTION 

In this section we derive the spectral function to first 
order in Poisson brackets, in the presence of a weak, time- 
dependent transverse driving field h(T). Since the spec- 
tral functions appearing in Eq. (jl6|) are already propor- 
tional to scattering rates, we will not keep gradient terms 
involving self energies in the following derivation. 

To proceed, we consider , which is determined from 
the relation 



where [G^]^^ = Gq ^ — E^, and G^^ is the inverse Green 
function in the clean limit, viz. Eq. (jlOp . To first order 
in Poisson brackets, one can show that in the Wigner 
representation 



where we once more have used the convenient matrix 
notation introduced in Eq. (jl5p . 

A similar relation can also be found for G'^, and we 
finally use that A = i(G^ - G^) to obtain an expression 
for the spectral function linear in gradients, viz. 



A{T;k,e)^Ao{T;k,s) 



where 



is the spectral function to zeroth order in gradients. Us- 
ing that we can rewrite 



we find the spectral function 



G«(r;fc,e)= l^i-- 

where now 

G^(T;fe,£) = 
1 



Gn AG 



G? + 0([...©, 



det (G^)-l V |7?>'i++S? £-e,-iA-Sfi, 



A(T;/c,e) = Ao(r;fc,£) 

+ ^:i^iA,(k,e)-r^dMk,E)\a,dT(h-a). (Al) 



Hence, the weak, transverse driving field induces only 
off-diagonal gradient contributions to the spectral func- 
tion. These gradient corrections prove to be essential in 
order to correctly capture the transverse magnetization 
dynamics. 
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